AFRL-RX-WP-TP-2008-4339 


NUMERICAL  AND  ANALYTICAL  MODELING  OF  LASER 
DEPOSITION  WITH  PREHEATING  (PREPRINT) 

Zhiqiang  Fan,  Jacquelyn  K.  Stroble,  Jianzhong  Ruan,  Todd  E.  Sparks,  and  Frank  Liou 
University  of  Missouri-Rolla 


MARCH  2007 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
MATERIALS  AND  MANUFACTURING  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7750 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 

Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY) 

March  2007 

2.  REPORT  TYPE 

Conference  Paper  Preprint 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

NUMERICAL  AND  ANALYTICAL  MODELING  OF  LASER  DEPOSITION  WITH 
PREHEATING  (PREPRINT) 

5a.  CONTRACT  NUMBER 

FA8650-04-C-5704 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

7801  IF 

6.  AUTHOR(S) 

Zhiqiang  Fan,  Jacquelyn  K.  Stroble,  Jianzhong  Ruan,  Todd  E.  Sparks,  and  Frank  Liou 

5d.  PROJECT  NUMBER 

2865 

5e.  TASK  NUMBER 

25 

5f.  WORK  UNIT  NUMBER 

25100000 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Missouri- Rolla 

1870  Miner  Circle 

Rolla,  MO  65409-0970 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

Materials  and  Manufacturing  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7750 

Air  Force  Materiel  Command 

United  States  Air  Force 

10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/RXLMP 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-RX-WP-TP -2008-4339 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

Conference  paper  submitted  to  the  proceedings  of  the  2007  ASME  International  Manufacturing  Science  and  Engineering 
Conference  (MSEC) ,  Atlanta,  GA,  15-17  Oct  2007.  PAO  Case  Number:  MLLM-0025;  Clearance  Date:  22  Mar  2007. 
This  work  was  funded  in  whole  or  in  part  by  Department  of  the  Air  Force  contract  FA8650-04-C-5704.  The  U.S. 
Government  has  for  itself  and  others  acting  on  its  behalf  an  unlimited,  paid-up,  nonexclusive,  irrevocable  worldwide 
license  to  use,  modify,  reproduce,  release,  perform,  display,  or  disclose  the  work  by  or  on  behalf  of  the  U.S. 

Government.  Paper  contains  color. 

14.  ABSTRACT 

Laser  deposition  allows  quick  fabrication  of  fully-dense  metallic  components  directly  from  CAD  solid  models.  This 
work  uses  both  numerical  and  analytical  approaches  to  model  the  laser  deposition  process  including  actual  deposition 
and  preheating.  The  numerical  approach  is  used  to  simulate  the  coupled,  interactive  transport  phenomena  during  actual 
deposition.  The  numerical  simulation  involves  laser  material  interaction,  free  surface  evolution,  and  melt-pool  dynamics. 
The  analytical  approach  is  used  to  model  heat  transfer  during  preheating.  The  combination  of  these  two  approaches  can 
increase  computational  efficiency  with  most  of  the  phenomena  associated  with  laser  deposition  modeled.  The  simulation 
is  applied  to  Ti-6A1-4V  and  simulation  results  are  compared  with  experimental  results. 


15.  SUBJECT  TERMS 

laser  deposition,  CAD,  numerical  modeling,  analytical  modeling,  preheating,  Ti-6A1-4V 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT: 

SAR 

18.  NUMBER 

OF  PAGES 

24 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

Todd  J.  Turner 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


i 


Numerical  and  Analytical  Modeling  of  Laser  Deposition  with  Preheating 

Zhiqiang  Fan,  Jacquelyn  K.  Stroble,  Jianzhong  Ruan,  Todd  E.  Sparks,  and  Frank  Liou 
Department  of  Mechanical  and  Aerospace  Engineering,  University  Of  Missouri-Rolla 
1870  Miner  Circle,  Rolla,  MO  65409 


Abstract 

Laser  deposition  allows  quick  fabrication  of  folly-dense  metallic  components  directly  from  CAD  solid 
models.  This  work  uses  both  numerical  and  analytical  approaches  to  model  the  laser  deposition 
process  including  actual  deposition  and  preheating.  The  numerical  approach  is  used  to  simulate  the 
coupled,  interactive  transport  phenomena  during  actual  deposition.  The  numerical  simulation  involves 
laser  material  interaction,  free  surface  evolution,  and  me  It -pool  dynamics.  The  analytical  approach  is 
used  to  model  heat  transfer  during  preheating.  The  combination  of  these  two  approaches  can  increase 
computational  efficiency  with  most  of  the  phenomena  associated  with  laser  deposition  modeled.  The 
simulation  is  applied  to  Ti-6A1-4V  and  simulation  results  are  compared  with  experimental  results. 


1.  Introduction 

Laser  deposition  is  an  extension  of  the  laser  cladding  process.  This  laser  additive  manufacturing 
technique  allows  quick  fabrication  of  folly-dense  metallic  components  directly  from  Computer  Aided 
Design  (CAD)  solid  models.  The  applications  of  laser  deposition  include  rapid  prototyping,  rapid 
tooling  and  part  refurbishment.  As  shown  in  Fig.  1,  laser  deposition  uses  a  focused  laser  beam  as  a 
heat  source  to  create  a  melt  pool  on  an  underlying  substrate.  Powder  material  is  then  injected  into  the 
melt  pool  through  nozzles.  The  incoming  powder  is  metallurgically  bonded  with  the  substrate  upon 
solidification.  The  part  is  fabricated  in  a  layer  by  layer  manner  in  a  shape  that  is  dictated  by  the  CAD 
solid  model. 
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Fig.  1.  Schematic  of  a  typical  laser  deposition  system 


During  the  laser  deposition  process,  several  defects,  such  as  porosity  and  cracks,  should  be  paid 
attention  to.  Cracks  initiate  corrosion  fracture  and  reduce  fatigue  strength  of  the  deposited  parts. 
Cracks  are  caused  by  the  residual  stresses  created  by  the  high  thennal  gradient  built  up  during  the 
cooling  stage.  Residual  stresses  can  be  reduced  by  a  reduction  of  the  cooling  rate.  This  can  be 
achieved  by  preheating  the  substrate.  Moreover,  the  preheating  causes  a  better  absorption  of  the  laser 
beam  and  further  it  is  possible  to  melt  more  powder  in  the  larger  melting  pool  and  enhance  the 
bonding.  Usually  the  preheating  is  performed  in  such  a  manner  that  a  very  small  melting  of  the 
substrate  material  occurs  without  powder  injection. 


Thennal  analysis  of  the  laser  deposition  process  is  very  important  for  optimization  of  the  process. 

If  the  substrate  surface  temperatures  remain  too  low,  wetted  deposit  material  is  limited.  In  that  case, 
irregularly  shaped  tracks  with  a  lot  of  cracks,  porosity  and  a  poor  bonding,  are  produced.  However,  if 
too  high  substrate  surface  temperatures  are  reached,  severe  melting  of  the  substrate  occurs.  The  high 
degree  of  dilution  can  deteriorate  the  clad  properties.  A  comprehensive  numerical  model  has  been 
developed  that  allows  the  prediction  of  temperature  distribution  and  melt  pool  dynamics.  This  model 
simulate  the  coaxial  laser  deposition  process  with  powder  injection,  and  considers  most  of  the 
associated  phenomena,  such  as  melting,  solidification,  evaporation,  evolution  of  the  free  surface,  and 
powder  injection.  Input  parameters  for  this  model  are  laser  machining  parameters  and  properties  of  the 
laser  beam,  as  well  as  material  properties  and  the  laser  beam  absorption.  To  get  more  accurate 
predictions,  finer  grid  needs  to  be  used.  This,  together  with  the  iterative  nature  of  the  numerical 
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algorithms,  causes  the  model  computationally  not  very  efficient.  In  this  paper  an  analytical  model  is 
applied  to  the  preheating  process  (without  powder  injection)  to  increase  the  computational  efficiency, 
while  the  actual  deposition  process  with  powder  injection  still  uses  the  numerical  model.  The  outputs, 
the  temperature  distribution  of  the  substrate,  are  used  as  the  initial  conditions  of  the  numerical  model. 

In  this  study,  a  coaxial  diode  laser  deposition  system,  LAMP  (developed  at  UMR),  is  considered  for 
simulations  and  experiments.  The  blown  powder  method  is  used  to  deliver  powder.  As  the  primary 
laser  in  the  LAMP  system,  the  diode  laser  is  believed  to  have  a  number  of  process  advantages  as 
opposed  to  the  CCL  or  Nd:YAG  laser,  including  the  better  material  coupling  efficiency  (laser 
absorption)  and  better  beam  profile  (non  Gauss-shape  intensity  distribution)  for  laser  deposition. 
Material  of  both  powder  and  substrates  is  Ti-6A1-4V,  which  is  widely  used  in  the  aerospace  industry. 

The  rest  of  the  paper  is  organized  as  follows:  In  Section  2,  the  numerical  model  is  presented.  In 
Section  3,  the  analytical  model  is  presented.  In  Section  4,  simulation  results  are  presented  and 
compared  to  the  experimental  results.  Discussions  are  made  in  this  section.  In  Section  5  conclusions 
are  completed. 

2.  Numerical  Model 

2.1  Governing  equations 

Figure  2  shows  a  schematic  diagram  of  the  calculation  domain,  including  the  substrate,  melt  pool, 
remelted  zone,  deposited  layer  and  part  of  the  gas  region.  In  the  laser  deposition  process,  melting  and 
solidification  cause  the  phase  transformation  at  the  solid/liquid  interface.  A  mushy  zone  containing 
solid  and  liquid  is  formed.  In  this  study  the  continuum  model  [1,2]  is  adopted  to  derive  the  governing 
equations. 
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The  assumptions  for  the  system  of  governing  equations  include:  (1)  the  fluid  flow  in  the  melt 
pool  is  a  Newtonian,  incompressible,  laminar  flow;  (2)  the  solid  and  liquid  phases  in  the  mushy  zone 
are  in  local  thermal  equilibrium;  (3)  the  solid  phase  is  rigid;  and  (4)  isotropic  permeability  exists.  For 
the  system  of  interest,  the  conservation  equations  for  mass,  momentum  and  energy  are  summarized  as 
follows: 

Continuity 

d^  +  V-(pV)  =  0  (!) 

dt 

Momentum 

-l(pf)  +  V-(pm  =  V-(M—Vn-Vp-^^(f-fs)  +  pg  +  S  (2) 

dt  p,  K  p, 

Energy 

+  V  •  (pf  h)  =  V  •  ( kVT )  -  V  •  (p(h,  -  h)(V  -  Vs )  (3) 

dt 

In  equations  (l)-(3),  the  continuum  density,  specific  heat,  thermal  conductivity,  vector  velocity,  and 
enthalpy  are  defined  as  follows: 

P  =  SsPs  +  S,Pi  c  =  fscs  +  f,c,  k  =  gsks  +  glkl 

V  =  ffs+fV,  h  =  fshs+fA  (4) 

The  liquid  fraction  temperature  relationship  is  given  by: 


2.  Schematic  of  the  calculation  domain 
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(5) 


0  if  T  <  Ts 
T-T 

- '  ifr<r<7; 

Tj-Ts 

1  if  T  >  T, 

The  other  volume  and  mass  fractions  can  be  obtained  by: 


f,= 


_S,P, 


f.= 


SsPs 


ss+g,=  1 


(6) 


P  P 

The  phase  enthalpy  for  the  solid  and  the  liquid  can  be  expressed  as: 

h,  =[c,(T)dT  h,  =  \T' c,(T)dT  +  f  c,(T)dT  +  L„ 


(7) 


where  Lm  is  the  latent  heat  of  melting. 

Permeability,  K,  is  assumed  to  vary  with  liquid  volume  fraction  according  to  the  Carman 
Kozeny  equation[3]  derived  from  Darcy’s  law: 


K  = 


S, 


c(i-g,y 


(8) 


where  the  parameter  C  is  a  constant  depending  on  the  morphology  and  size  of  the  dendrites  in  the 
mushy  zone.  The  S  is  a  source  term  that  will  be  defined  in  Section  2.4. 


2.2  Tracking  of  the  solid/liquid  interface 


The  solid/liquid  interface  is  implicitly  tracked  by  the  continuum  model  [1,  2],  In  the  solid  phase 
region  and  liquid  phase  region,  the  third  term  on  the  right-hand  side  of  Eq.  (2)  vanishes.  This  is 
because  in  the  solid  phase  region  V  =  Vs  =  0  and  in  liquid  phase  region  K— »<x>  since  gi  =  1 .  So  this 
term  is  only  valid  in  the  mushy  zone. 

2.3  Tracking  of  the  free  surface 

The  liquid/vapor  interface,  or  the  free  surface  of  the  melt  pool,  is  very  complex  due  to  surface 
tension,  thermocapillary  force,  and  impaction  of  the  powder  injection.  In  this  study,  the 
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Volume-Of-Fluid  (VOF)  method  [4]  is  employed  to  track  the  evolution  of  the  moving  free  surface  of 
the  melt  pool.  The  melt  pool  configuration  is  defined  in  terms  of  a  volume  of  fluid  function,  F(x,y,t), 
which  represents  the  volume  of  fluid  per  unit  volume  and  satisfies  the  conservation  equation: 

Z7 

—  +  (V-V)F  =  0  (9) 

dt 


2.4  Formulation  of  source  term 


The  source  term,  S,  in  the  momentum  equation  is  contributed  by  the  interface  forces  acting  on  the 
free  surface,  such  as  surface  tension,  etc.  In  this  study,  the  continuum  surface  force  (CSF)  model  [5]  is 
used  to  refonnulate  the  surface  force.  In  its  standard  form,  surface  tension  is  fonnulated  as  [5]: 


Fs{xs)^hyK  +  V  sy 


(10) 


where  Fs  (xs )  is  the  net  surface  force  at  a  point  xs  on  a  interface  S.  h  is  a  unit  normal  to  S  at  the 


point  x ,  which  is  given  by: 


n 

n~V~\ 

\n\ 


(ID 


where  Tl  is  the  surface  normal  vector  and  can  be  computed  from  the  gradient  of  the  VOF  function: 

fl  =  VF  (12) 

V5  is  the  gradient  along  a  direction  tangential  to  the  interface,  which  is  defined  as: 

Vs  =V-V7v  =V- n(n  -V)  (13) 

y  and  k  represent  surface  tension  coefficient  and  curvature,  respectively,  k  is  given  in  [6]: 


K  =  -(V  •  n)  =  — 
n 


(r-v)|»|-(v.») 


(14) 


By  using  the  CSF  model,  the  surface  force  Fs  is  reformulated  into  a  volume  force  Fh  as  follows: 


Fb  =  Fs  n 


F 

F 


clVFl 


F_ 

F 


(15) 
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where  ( F \  is  the  averaged  F  value  across  the  free  surface.  Thus,  the  source  tern  S  in  equation  (2) 


formulated  as: 


S  =  (n/K  +  V5^)|VF| 


2.5  Boundary  Conditions 


The  boundary  conditions  at  the  free  surface  satisfie  the  following  equation: 


ST  _  Glaser  -Pauen) 


~K(T  ~  Tx)  -  sg(Ta  —T*)  —  meLv 


where  terms  on  the  right-hand  side  are  laser  irradiation,  convective  heat  loss,  radiation  heat  loss  and 
evaporation  heat  loss,  respectively.  PiaSeris  the  power  of  laser  beam,  Patten  the  power  attenuated  by  the 
powder  cloud,  R  is  the  laser  beam  radius,  q  the  laser  absorption  coefficient.  Patten  is  calculated 
according  to  Frenk’s  et  al.  model  [7]  with  a  minor  modification: 


Patten  =  Plaser  1  “  eXP  “ 


}QeJ± 

7tprD.vp 


where  til  denotes  the  powder  mass  flow  rate,  1  is  the  stand-off  distance  from  the  nozzle  exit  to  the 
substrate,  p  is  powder  density,  rp  is  the  radius  of  the  powder  particle,  Djet  is  the  diameter  of  the  powder 
jet,  vp  is  the  powder  injection  velocity,  and  Qext  is  the  extinction  coefficient,  ft  is  assumed  that  the 
extinction  cross  section  is  close  to  the  actual  geometrical  cross  section,  and  Qext  takes  a  value  of  unity. 

In  the  evaporation  term,  tii  is  the  evaporation  mass  flux  and  Lv  is  the  latent  heat  of  evaporation. 


According  to  Choi  et  al.’s  “overall  evaporation  model”  [8],  fh  is  of  the  fonn: 


1  8876 

\ogme  =A  +  6.1210 - ^-0.51ogr 


where  A  is  a  constant  dependent  on  the  material. 


The  boundary  conditions  at  the  bottom,  left  and  right  wall  satisfy  the  following  equations: 
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(21) 


tST 

8n 


-K(T~TJ 


u  =  0  v  =  0 


(22) 


2.6  Powder  injection 

Powder  particles  that  inject  onto  the  top  surface  can  be  classified  into  three  categories.  (1)  Those 
powder  particles  that  have  not  been  melted  during  their  passage  and  hit  the  solid  part  of  the  substrate 
will  deflect  and  lose.  (2)  Those  powder  particles  that  have  been  melted  before  they  arrive  on  the 
substrate  and  impact  the  solid  part  of  the  substrate  will  stick  to  the  surface  of  the  substrate.  (3)  Those 
powder  particles  that  fall  into  the  melt  pool  including  melted  and  unmelted  will  merge  and  mix  with 
the  molten  liquid  in  the  melt  pool  [9].  In  this  study,  the  powder  particles  belonging  to  the  last  two 
cases  will  be  utilized. 

2.7  Numerical  algorithm 

The  governing  equations  (Eq.  (l)-(3),  and  (10))  and  all  related  supplemental  equations  and  boundary 
conditions  are  solved  through  the  following  iterative  scheme: 

1.  Equation  ( l)-(2)  and  the  related  boundary  conditions  are  solved  iteratively  using  the  SOLA-VOF 
algorithm  [10]  to  obtain  velocities  and  pressures.  At  each  time  step,  the  discretized  momentum 
equations  calculate  new  velocities  in  tenns  of  an  estimated  pressure  field.  Then  the  pressure  field 
is  iteratively  adjusted  and  velocity  changes  induced  by  each  pressure  correction  are  added  to  the 
previous  velocities.  This  iterative  process  is  repeated  until  the  continuity  equation  is  satisfied 
under  an  imposed  tolerance  by  the  newly  computed  velocities. 

2.  The  energy  equation  (4)  is  solved  by  an  implicit  method. 

3.  Equation  (10)  is  solved  to  obtain  the  updated  free  surface  and  geometry  of  the  melt  pool. 

4.  Advance  to  the  next  time  step  and  back  to  step  1  until  the  desired  process  time  is  reached. 

In  this  scheme,  staggered  grids  are  employed  where  the  temperatures,  pressures  and  VOF  function  are 
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located  at  the  cell  center  and  the  velocities  at  the  walls.  The  source  tenn  in  the  momentum  is  obtained 


at  grid  points  that  are  located  in  the  transition  region.  A  fixed  grid  is  used  with  the  mesh  size  of  10pm, 
which  is  proved  to  achieve  grid  independence.  The  time  step  is  taken  at  the  level  of  10"6s  initially  and 
adapted  subsequently  according  to  the  convergence  and  stability  requirements  of  the  SOLA-VOF 
algorithm  [10],  the  CSF  model  [5],  and  the  numerical  solution  for  the  energy  equation. 


3.  Analytical  Model 

In  this  paper  the  analytical  solution  for  the  temperature  distribution  in  preheating  of  the  substrate  is 
based  on  Komanduri  and  Hou’s  work  [11],  which  is  based  on  moving  heat  source  theory  of  Jaeger 
[12]  and  Carslaw  and  Jaeger  [13]. 


3.1  Solution  of  a  disc  heat  source  with  a  uniform  intensity  distribution,  moving  on  the  surface  of 
a  semi-infinite  medium 

Follow  Komanduri  and  Hou’s  work  [11],  the  temperature  rise  of  a  point  M(x,  y,  z)  at  any  time  t  after 
the  initiation  of  the  heat  source  in  a  semi-infinite  medium  due  to  a  disc  heat  source  is  given  by: 

P  (t  dz  ,  r2  +  X2  +  y2  +  z2  ,  .  Xv 


A T  = - : - —  f  -^-exp(- 

2cp(4n:a)  r 


)  •  exp(-A— )  •  exp(-^T)  •  I0(^(X  +vT)2+y2)  (23) 
4ar  2  a  4  a  2ar 


where  y,  z  are  coordinates  in  the  absolute  coordinate  system,  while  X  is  the  corresponding  moving 
coordinate.  P  is  the  laser  power  absorbed  by  the  substrate,  ro  is  the  radius  of  the  disc  heat  source,  v  its 
moving  velocity,  I0  the  modified  Bessel  function  of  the  first  kind,  order  zero,  c,  p,  a  are  specific  heat, 

density  and  thermal  diffusivity,  respectively.  Let  R  =  r0  +  X  +  y  +  z  ,  Eq.  (23)  becomes: 


A  T  = 


f  4tt exp(- )  •  exp(- Af )  •  exp(- •  70 yj(X  +vr)2  +  y2 ) 

JO  T  dr/r  In  dr/  1  nr 


2cp(A*7ia)  Jo  T 


jr 

4ar' 


Xv 
2  a 


(24) 


4  a ' 


2ar 


3.2  Solution  of  a  disc  heat  source  with  a  uniform  intensity  distribution,  considering  the 
boundary  conditions  of  the  bottom  and  side  surfaces 

Here  the  bottom  and  side  surfaces  are  considered  as  adiabatic.  Using  the  image  method,  nine  image 
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heat  sources  are  considered.  The  temperature  rise  at  any  point  M  is  the  sum  of  the  effects  from  all  the 


primary  and  the  image  heat  sources  which  are  located  away  from  the  point  M  at  distances  of  Ro,  Ri, 
R2,  R3,  R4,  R5,  Re,  R7,  Rs,  R9,  respectively.  Suppose  the  substrate  has  the  dimensions  H  x  L  x  W,  H 
is  the  height  of  the  substrate,  L  is  the  length  of  substrate  in  the  heat  source  moving  direction.  The 
values  of  these  distances  and  their  projections  on  the  X-axis  are  given  in  Table  1. 


Table  1  Distances  and  their  projections 


Symbol  of 

Value  of  Distance 

Symbol  of 

Value  of  Projection 

Ro 

V-^+r+z2 

X0 

X 

Ri 

7 X 2  +  V2  +  (2 H  -  z)2 

Xi 

X 

r2 

7(X  +  2  vt)2+y2  +z2 

x2 

-  (X  +  2vt) 

R3 

■J(X  +  2vt)2  +  y2  +  (2H-z)2 

x3 

-  (X  +  2vt) 

R4 

7(2  (L-vt)-X)2+y2+z2 

x4 

1 

> 

1 

Rs 

7 (2(L  -  vt)  -  X)2  +  v2  +  (2H  -  z)2 

X5 

2(L-vt)-X 

Re 

7  X2+(y  +  W)2+z2 

X6 

X 

Rt 

■yjx2  +  (v  +  W)2  +  (2H  -  z)2 

X7 

X 

Rs 

7  X2  +  (y  +  W)2+z2 

X8 

X 

r9 

■\Jx2  +  (v  +  W)2  +  (2H  -  z)2 

X9 

X 

The  temperature  rise  at  any  point  M  caused  by  the  primary  heat  source  is  given  by  Eq.  (24),  and  the 
temperature  rise  at  any  point  M  caused  by  each  of  the  image  heat  source  is  obtained  by  substituting 
the  aforementioned  values  of  the  distances  and  their  relevant  projections  on  the  X-axis  instead  of  R 
and  X  in  Eq.  (24). 


4.  Numerical  Simulation  and  Experiments 

4.1  Numerical  simulation 

Simulation  is  perfonned  based  on  the  capability  of  our  experimental  facilities  to  compare  the 
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simulation  results  with  the  experimental  measurements.  A  continuous  wave  diode  laser  with  an  808 
nm  wavelength  is  considered  as  the  energy  source.  The  laser  intensity  distribution  is  unifonn.  For 
substrates,  Ti-6A1-V4  plates  with  a  thickness  of  0.25  inch  are  selected.  Ti-6A1-V4  powder  particles 
with  a  diameter  from  40  to  140  pm  are  used  as  deposit  material.  The  laser  absorption  coefficient  is 
measured  by  Sparks  et  al.  [14].  The  material  properties  and  the  main  process  parameters  are  shown  in 
Table  2.  Figure  2  shows  the  typical  simulation  results  for  temperature,  velocity  and  VOF  function. 


(a)  Temperature  field  (b)  Velocity  field  (c)  Volume  of  Fluid 

Fig.  2.  Simulation  results  at  t  =  130  ms  (laser  power-9 10W,  travel  speed-  20ipm,  powder  mass  flow 
rate-4. 68g/min) 

4.2  Experiments 

The  experiments  were  performed  on  the  LAMP  system  shown  in  Fig.  3.  The  system  consists  of  a 
diode  laser,  powder  delivery  unit,  5-axis  CNC  machine,  and  monitoring  subsystem.  The  laser  system 
used  in  the  study  was  Nuvonyx  (Nuvonyx  Inc.)  ISL-1000M  Laser  Diode  System  which  combines 
state-of-the-art  micro-optics  with  laser  diodes  to  produce  the  only  single  wavelength  fiber  coupled 
direct  diode  laser  at  power  levels  up  to  1000  watts  CW.  The  laser  emits  at  808  nm  and  operates  in  the 
continuous  wave  (CW)  mode.  To  protect  oxidization  of  Ti-6A1-V4,  the  system  is  covered  in  an 
environmental  chamber  to  supply  argon  gas  for  titanium  deposition.  For  the  other  aspects  of  the 
system  architecture,  refer  to  Liou  et  al.  [15]  and  Boddu  et  al.  [16].  The  substrates  have  dimensions  of 
2.5x2. 5><0. 4  in.  The  Ti-6A1-4V  samples  were  irradiated  using  a  laser  beam  with  a  beam  spot  diameter 
of  2.5  mm. 
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Table  1.  Material  properties  for  Ti-6A1-4V  and  main  process  parameters 


Nonmenclature 

Symbol 

Value  (unit) 

Melting  temperature 

T 

1  m 

1900.0K 

Liquidus  temperature 

T, 

1923. OK 

Solidus  temperature 

Ts 

1877.0K 

Evaporation  temperature 

Tv 

3533. OK 

Solid  specific  heat  at  constant 
pressure  [17] 

Cps 

1 

483.04  +  0.215T  T  <1268.6: 

J/kgK 

412. 7  +  0. 180  IT  1268  <  T  <  1923 

Liquid  specific  heat  at  constant 

Cpl 

831.0  J/kgK 

pressure  [18] 

Thermal  conductivity  [17] 

k 

1.2595 +  0.0157T  T  <  1268AT 

< 

3.5127 +  0.0127T  1268<T<1923  W/mK 

-12.752 +  0.024T  T>1923 

Solid  density  [18] 

ps 

4420-0.154  (T-298K) 

Liquid  density  [18] 

Pi 

3920-0.68  (T-  1923  K) 

Latent  heat  of  fusion  [18] 

Lm 

2.86  x  105  J/kg 

Latent  heat  of  evaporation 

Lv 

9.83  x  106  J/kg 

Dynamic  viscosity 

P 

3.25  x  10'3  N/m  s  (1923K)  3.03  x  10'3 

(1973K) 

2.66  x  10'3  (2073K)  2.36  x  10'3  (2173K) 

Radiation  emissivity  [19] 

8 

0.1536  +  1.8377  x  10 4  (T  -  300.0  K) 

Laser  absorption  coefficient 

n 

0.4 

[14] 

Powder  particle  diameter 

Dr 

40-140  pm 

Shielding  gas  pressure 

Pg 

5  psi 

Ambient  temperature 

Too 

300K 

Convective  coefficient 

hc 

10  W/m2K 

In  order  to  validate  the  model  predictions,  single  path  deposition  experiments  are  conducted.  The 
comparisons  between  model  predictions  and  experimental  results  are  conducted  in  terms  of  melt  pool 
peak  temperature.  The  melt  pool  peak  temperature  is  calibrated  through  the  dual-wavelength 
non-contact  temperature  sensor,  which  can  effectively  decrease  the  disturbance  from  the  powder  and 
other  dusts. 
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CMOS  Camera 

Laser  ontic  fiber 

Shielding  gas 
Powder  and  direct  gas 
Lens 

_ A 

Substrate 

Fig.  3.  Schematic  of  experimental  setup 

4.3  Comparisons  and  discussions 

Fig.  4  and  Fig.  5  show  the  comparisons  between  experimental  measurements  and  model  predictions. 
Fig.  4  shows  the  effects  of  laser  power  on  melt  pool  peak  temperature.  It  can  be  seen  that  an  increase 
in  the  laser  power  will  increase  the  melt  pool  temperature.  This  is  easy  to  understand.  As  the  laser 
power  increases,  more  power  is  available  for  melting  the  substrate.  Fig.  4  shows  the  effects  of  laser 
scanning  speed  on  the  melt  pool  peak  temperature.  An  increase  in  the  laser  scanning  speed  will 
decrease  the  melt  pool  peak  temperature.  This  is  because,  as  scanning  speed  decreases,  the  laser 
material  interaction  time  is  extended. 

From  Fig.  4-5,  we  can  see  that  the  general  trend  between  experimental  measurements  and  model 
predictions  is  consistent.  At  different  power  intensity  level,  there  is  a  different  error  from  10  K  (about 
0.5%)  to  121K  (about  5%).  It  can  be  seen  that  at  a  higher  power  intensity  level,  there  is  a  bigger  error 
between  measurements  and  predictions.  This  is  because  the  numerical  model  is  two-dimensional.  It 
doesn’t  consider  the  heat  and  mass  transfer  in  the  third  direction.  At  a  higher  power  intensity  level, 
heat  and  mass  transfer  in  the  third  dimension  is  more  significant. 

The  errors  between  experimental  measurements  and  model  predictions  are  analyzed  to  mainly  come 
from  the  following  aspects:  (l)The  two-dimensional  nature  of  the  numerical  model;  (2)The 
thermo-physical  properties  taken  for  the  analytical  model;  (3)  The  uncertainties  of  the  material 
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properties  and  the  appropriateness  of  the  sub-models  for  the  numerical  model;  (4)  Boundary 
conditions.  Adiabatic  boundary  conditions  are  assumed  in  the  analytical  model  and  the  numerical 
model  for  the  bottom  surface  and  side  surfaces.  Measurements  also  have  been  taken  to  achieve  such 
boundary  conditions  in  experiments.  But  it  is  hard  to  get  absolute  adiabatic  boundary  conditions. 


Laser  Power,  W 

Fig.  4  Comparisons  between  experimental  measurements  and  model 
predictions  at  a  constant  powder  mass  flow  rate  of  5g/min  and  a 
constant  laser  scanning  speed  of  20  inch/min.  (For  actual  deposition) 


Laser  Scanning  Speed,  inch/min. 

Fig.  5  Comparisons  between  experimental  measurements  and  model 
predictions  at  a  constant  powder  mass  flow  rate  of  5g/min  and  a 
constant  laser  power  of  700  W.  (For  actual  deposition) 
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5.  Conclusions 


An  analytical  model  for  thennal  analysis  of  temperature  rise  due  to  a  moving  heat  source  is  combined 
into  a  comprehensive  heat  transfer  and  fluid  flow  numerical  model  for  the  laser  deposition  process. 
The  analytical  model  is  used  for  the  preheating  process  before  the  actual  laser  deposition  with  powder 
injection.  And  the  numerical  model  is  used  for  the  actual  laser  deposition  process.  Thus  the  outputs  of 
the  analytical  model  are  used  as  the  inputs  of  the  numerical  model.  Experiments  have  been  done  to 
validate  the  model  predictions.  A  consistent  general  trend  is  found  between  experimental 
measurements  and  the  model  predictions.  The  sources  of  the  errors  have  been  analyzed. 
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